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Abstract 

Salmonella isolates harbour a range of resident prophages which can influence their virulence and ability to compete and 
survive in their environment. Phage gene profiling of a range of phage types of Salmonella enterica subsp. enterica serovar 
Typhimurium (S. Typhimurium) indicates a significant level of correlation of phage gene profile with phage type as well as 
correlation with genotypes determined by a combination of multi-locus variable-number tandem repeat (VNTR) typing and 
clustered regularly interspaced short palindromic repeats (CRISPR) typing. Variation in phage gene profiles appears to be 
partly linked to differences in composition of variants of known prophages. We therefore conducted a study of the 
distribution of variants of ST64B and Gifsy-1 prophages and coincidently the presence of Gifsy-3 prophage in a range of S. 
Typhimurium phage types and genotypes. We have discovered two variants of the DTI 04 variant of ST64B and at least two 
new variants of Gifsy-1 as well as variants of related phage genes. While there is definite correlation between phage type 
and the prophage profile based on ST64B and Gifsy-1 variants we find stronger correlation between the VNTR/CRISPR 
genotype and prophage profile. Further differentiation of some genotypes is obtained by addition of the distribution of 
Gifsy-3 and a sequence variant of the substituted SB26 gene from the DTI 04 variant of ST64B. To explain the correlation 
between genotype and prophage profile we propose that suites of resident prophages promote clonality possibly through 
superinfection exclusion systems. 
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Introduction 

Strains of Salmonella enterica serovar Typhimurium may contain a 
range of temperate phages residing as prophages in their 
chromosomes. The prophage load of the host bacterium can be 
very important to its ability to compete and survive in its 
environment because prophage genomes may contain virulence 
and fitness factors which become a driving force in the pathogen- 
host interaction and can lead to emergence of new epidemic clones 
[1]. Furthermore, low level release of bacteriophage from 
lysogenised strains may give them a competitive advantage over 
phage-sensitive strains [2] . The phage type of a strain is likely to 
reflect the prophage composition of the isolate because the lysis 
patterns of the phages in the typing panel may be altered by new 
phage acquisitions or deletions [3] . Knowledge of the prophage 
composition as well as the underlying genotype of the bacterial 
host can give clues to the possible phylogenetic origins of newly 
emerged epidemic strains [2]. 

According to current theory, bacteriophage genomes are a 
mosaic consisting of combinations of modules which are 
exchangeable among members of the bacteriophage population 
by either homologous or illegitimate recombination events [4]. A 
comparison of related phage genomes will typically show regions 



of close sequence identity interrupted by regions of non-identity. 
These gaps are mainly DNA replacements where the DNA 
segment in one phage is replaced in another phage by a segment of 
unrelated DNA that frequently fulfils the same or related function 
[!]• 

The ST64B bacteriophage (henceforth termed ST64Bdts4) 
induced from S. Typhimurium DT64 [3] is described as a genetic 
mosaic which has acquired significant portions of its genome from 
sources outside the genus Salmonella [5]. ST64B DX64 is grouped 
with P27-like phages [6] but has similarities to lamboid phages. 
Several studies have made limited comparisons of the ST64B 
prophages in different S. Typhimurium isolates. Figueroa-Bossi 
and Bossi [7] reported that the ST64B prophage in DTI 04 
showed some sequence differences in the SB21-22 region 
compared with ST64B from ATCC 14028s and SL1344. Ross 
and Heuzenroeder [8] reported that the SB26 gene was the least 
common of targeted ST64Box64 genes in a panel of isolates which 
included DTs 126, 108, 170, 12 and 12a but not DTI 04. Hermans 
et al [9] using primer sets derived from S. Typhimurium DTI 04 
strain 7945 showed that some S. Typhimurium strains tested 
positive for presence of ST64B but lacked the internal sequence 
from the DT 1 04-derived ST64B. Cooke et al [10] used 
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microarray analysis and PGR targeting prophage sequences to 
differentiate a panel of 23 S. Typhimurium isolates, mainly DTI 04 
and DTI 70. The primers for ST64B were specific for a sequence 
in the DTI 04 form of ST64B (henceforth termed ST64B DT104 ) 
which could not be amplified in ST64Bm- 64 but this was not 
recognized at the time. 

Figueroa-Bossi et al [1 1] first described the closely related Gifsy- 
1 and Gifsy-2 lamboid prophages present in both LT2 and ATCC 
14028 strains of S. Typhimurium and subsequently reported a 
third type of related lamboid phage in ATCC 14028 which they 
called Gifsy-3 [12]. Gifsy-1 and Gifsy-2 occur commonly in 
serovar Typhimurium but Gifsy-3 is found only rarely [13]. Whilst 
Gifsy-2 is highly conserved within serovar Typhimurium, Gifsy-1 
shows extensive strain-to-strain variability within serovar Typhi- 
murium [13,14]. DT104 NCTC 13348 carries a Gifsy- 1 -related 
prophage, Gifsy- 1dtio4, an d DTI 04 is not immune to infection by 
Gifsy-1 from ATCC 14028 [2]. Gifsy-1 from LT2 ATCC 43971 
(Gifsy-lLT2) contains the GipA moron [1] but Gifsy-lDTi04 lacks 
GipA [13]. Gifsy-1 from LT2, ATCC 14028 and SL1344 all have 
different immunity modules [15]. Gifsy- 1 L t2 nas the same 
immunity module as Gifsy-2, while Gifsy- 1 in SL1344 (Gifsy- 
1sli344) has the immunity of Gifsy-3 from ATCC 14028 (Gifsy- 

3 14028)- 

Phage typing of S. Typhimurium according to the scheme of 
Anderson [16] has been used for many years to further 
characterise isolates for epidemiological purposes. The phage 
sensitivity patterns (or phage type) are determined by the 
properties of the infecting phages, cell surface factors on the 
isolate and the susceptibility of the typing phage to immunity-type 
repressors and to various superinfection exclusion systems 
controlled by residing prophages. Phage type can be altered by 
the introduction into the bacterial cell of mobile elements 
particularly bacteriophages [3] and plasmids [17]. The relation- 
ships among different phage types are poorly understood. 

MLVA (multilocus variable-number-tandem-repeats analysis) 
typing is technically less demanding than phage typing and has 
been found to give generally superior discrimination between 
isolates [18]. It is more informative than phage typing because it 
allows isolates to be compared genetically at the chromosomal 
level and possible phylogenetic relationships to be inferred. The 
Public Health Microbiology Laboratory at Forensic and Scientific 
Services has used MLVA for typing of all isolates of S. 
Typhimurium since 2006 and has more recently supplemented 
this with additional tandem repeats and with clustered regularly 
interspaced short palindromic repeats (CRISPR) analysis (http:// 
crispr.u-psud.fr/crispr/). The combination of two typing methods 
has allowed recognition of distinct genotypes termed Repeats 
Groups (RGs) which correlate strongly with single or numbers of 
phage types. Trials with phage gene profiling of S. Typhimurium 
isolates have also been conducted [19]. These have included PCR 
tests for five genes found in the phage ST64B DT64 , namely SB6, 
SB26, SB28, SB37 and SB46 (CC-5). We noted that some 
genotypes including the one which included DTI 04 NCTC 13348 
lacked SB26 but had the other ST64B genes while other genotypes 
had only the SB46 gene. We hypothesised that the isolates missing 
only SB26 gene had the ST64B D xio4 prophage. Investigation into 
why some genotypes had only SB46 gene led to the discovery of a 
second set of ST64B genes in Gifsy-lDTio4 but not in Gifsy-lLT2- 
An oblique reference to this phenomenon was later noted in 
Figure 3 from Lemire et al [13]. 

There is a need for a more thorough investigation into the 
existence of variants of ST64B and Gifsy-1 prophages. We 
conducted a more comprehensive comparison of available 
prophage sequences in order to find specific sequences which 



could identify recognized variants and carried out a survey of well 
characterised isolates to see how variants are distributed amongst 
the different genotypes of S. Typhimurium. 

Results 

Comparison of ST64B Sequences in DT64, DT2 and 
DT104 

The ST64B sequences from DT64 (ST64B DT64 ), DT2 
(ST64B DT2 ) and DT104 (ST64B DT104 ) were 40149, 40147 and 
39463 bp respectively in total length. The overall differences 
between ST64B DT64 and ST64B DT2 sequences were confined to 
relatively few single-nucleotide polymorphisms (SNPs) and Indels 
so that comparison of ST64BrjTi04 with ST64Bdt2 was essentially 
the same as with ST64BrjT64- The comparison of ST64Bdtio4 an d 
ST64B DX2 showed that there were six regions of non-identity 
corresponding to ST64B DT64 genes SB32-34, SB37-38, SB49, 
SB51, SB53-55, SB26+half SB27 and one of only 66% identity 
(SB52) (Figure S3). A visual comparison of the two ST64B 
sequences using MAUVE is shown (Figure 1). 

Testing of Isolates for ST64B DT64 and ST64B DT104 
Sequences 

Distribution of SB26 gene of ST64B DT6 4. We tested for 
SB26 gene of ST64Bdt64 m all 214 isolates belonging to the panel 
of well characterised isolates (Figures SI and S2). These have been 
divided into Major Groups 1 and 2 with each Major Group 
subdivided into numbers of Repeats Groups. We found that SB26 
was detected in all but one of the isolates in Major Group 2 but 
only in two RDNC isolates belonging to Repeats Group 2A 
(RG2A) in Major Group 1. The DNA sequences of the SB26 allele 
in these two RDNC strains and for a few controls from isolates in 
Major Group 2 were the same as for the SB26 allele in 
ST64BrjT64- On the basis of the phage gene profiling results 
[19] and testing for the SB26 gene in the characterised panel of 
isolates we made the assumption that isolates positive for the SB26 
gene probably contain the ST64B DT64 prophage. 

RG2 (2A+2B) consists of mostly DT4, DT141 and DT141 var 1 
and 2 isolates and includes the sequenced strain LT2 ATCC 
43971. All members of RG2 have an STTR11 allele of 716 bp 
rarely seen in isolates of any other Repeats Group. We tested a 
further 59 isolates which had MLVA types consistent with those of 
RG2 and the 716 bp STTR11 allele. Eleven isolates tested 
positive for the SB26 gene and were positive also for SB6, SB28, 
SB 3 7 and SB46 genes of ST64B DT64 . The presence of all five 
genes was taken as presumptive evidence that the whole 
ST64B DX64 prophage was present in these isolates. Further 
examination showed that all the RG2 isolates with the SB26 gene 
belonged to RG2A which is distinguished from RG2B on the basis 
of different CRISPR 2 profiles (Figure S2). 

Distribution of ST64Bi>Ti04-specific sequences. We vali- 
dated the PCR tests for the ST64B DT104 -specific sequences by 
testing isolates from RG1 (DT170 and DT12) and RG8 (DT104L 
and DT 1 2L) which were expected to be positive from phage gene 
profile results [19] and an isolate from RG13A expected to be 
negative. RG1 and RGS isolates were positive for all seven 
ST64B D xio4-specific sequences and the RG13A isolate was 
negative for all seven sequences. All seven PCR products for one 
isolate from each of RG1 and RG8 were sequenced and found to 
have 100% identity with the DTI 04 sequences. This was 
presumptive evidence that the samples with all seven 
ST64B D Tio4-specific sequences harboured the whole ST64B DT104 
prophage. 
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Figure 1. Comparison of ST64B sequences for 5. Typhimurium strains DT2 (Sanger) and DT104 (NCTC 13348) by MAUVE (See also 
Figure S3 for gene locations and levels of identity for various regions). 

doi:1 0.1 371 /journal.pone.0086203.g001 



We subsequently undertook testing of a representative selection 
of isolates from the various RGs in Major Groups 1 and 2 for all 
seven ST64B DT io4 phage sequences. None of the Major Group 2 
isolates had any of the ST64Bdtio4 sequences. All of the isolates in 
RG1 and RGS were positive for all of the ST64B DT104 -specific 
sequences. One of four isolates in RG6A with a 171-X-0-0-462 
MLVA profile type was positive for the seven sequences and all 
four of the isolates in RG6B with a 1 7 1 -x-x-0-489 MLVA type 
were positive as well. An isolate typed anomalously as DTI 26 but 
with MLVA and CRISPR types consistent with membership of 
RG2A was positive for all seven ST64B D tio4 sequences. All 
isolates in RG9A were positive only for the SB26dtio4 gene 
sequence but isolates in RG9B were negative for this sequence. 
RG9A consists mostly of DTI 97 isolates with a small number 
typed as DT6 and DT43, and one, anomalously, as DT104L 
whereas RG9B consists of DTI 2a isolates. RG9B is distinguished 
from RG9A by a significantly different CRISPR 2 profile (Figure 
S2) although the CRISPR 1 profile is identical to that of RG9A 
(Figure SI). The SB26dtio4 gene product from RG9A isolates was 
found to have a nucleotide sequence which shared only 97% 
identity with the sequence in ST64B DT io4- In RG7, two isolates, 
05P41568972 typed as DTI and 09P84391818 as DT41, with the 
same MLVA and CRISPR profiles, along with ATCC 13311, 
which had almost the same MLVA profile, were all found to be 
positive for the ST64B D xi04 sequences corresponding to SB26 
gene, SB33 gene and SB37-38 intergenic sequence and for SB28 
and SB 3 7 shared by both types of ST64B prophage but not for the 
ST64B DTI( )4-specific SB49, 50/51 intergenic, SB52 and SB53 
sequences and not for SB6 or SB46 common to both ST64B 
prophages. It appears therefore that these isolates have a variant of 
ST64B DTH )4 in which the sequence between SB37 to beyond SB6 
is altered to some extent. 

We wished to clarify the ST64B status of RG6B members with a 
171-x-x-0-489 MLVA profile which is shared by members of RG3 
(Figures SI and S2) by testing a larger number of isolates 
belonging to the two groups. Thirty isolates assigned to RG3 
consisting of mostly DT101, DTI 2a and a few RDNC isolates 
with an STTR7 allele of 620 bp (except for two with 582 bp allele) 
and an ST3 allele of 200 bp (except for two with 189 bp allele) 
were found to be uniformly negative for ST64B D xi04-specific 
sequences. Forty isolates assigned to RG6B consisting of DTI 20 
(8), DT193 (12), DT195 (1), U302 (2), RDNC (2) and Untypable 
(15) isolates, which aU had STTR7 alleles of 582 bp and ST3 



alleles of 189 bp, were all positive for ST64B D xio4-specific 
sequences. It will be shown later in the results section, that the 
RG3 isolates have a different Gifsy-1 prophage from that in 
RG6B. It was also noted from the antibiotic sensitivity results that 
accompany phage typing results that 29 out the 30 isolates in RG3 
had a fully sensitive antibiotic profile whereas 38 of the 40 isolates 
in RG6B had a multi-resistant antibiotic profile most commonly 
AMP STR TET SUL. Significantly too, clinical notes supplied 
with the request slips indicated that at least 8 patients with RG6B 
isolates had a history of travel to SE Asia (7) or Middle East (1) 
while there was no travel history for any patients with the RG3 
isolates. Furthermore, RG6B isolates had not been seen prior to 
2006 and the numbers have been increasing every year since then 
whereas the RG3 isolates had been much more common before 
2006 and numbers had fallen to very low levels since that time. As 
a consequence of the testing for ST64B D xio4 we have been able to 
recognise these isolates as an emerging S. Typhimurium genotype 
possibly of Asian origin. 

To further explore the status of RG9A members we selected 
from the years 2008 and 2009 thirty seven DT197 isolates with a 
variety of MLVA profiles for amplification of two contiguous 
DT 1 04-specific SB26 sequences of 585 and 897 bp. All 37 isolates 
were positive for both sequences. Sequencing of both products for 
three isolates with different MLVA profiles showed that all three 
had identical sequences. The combined 1482 bp allele had 26 
SNPs, a single nucleotide insertion and a single nucleotide deletion 
relative to the SB26 sequence in ST64Brjxio4 (GenBank Accession 
No. KC 172925.1). The DTI 04-specific SB26 status of earlier 
DTI 97 isolates was also investigated (Text SI). 

Variation in SB46 gene sequences. Testing for the SB46 
gene common to both ST64B DT 64 and ST64B DT10 4 showed that 
all selected members of RG3, RG4 and RG9 comprising mainly 
DTI 2a, DT101, DT179 and DT197 isolates were positive for 
SB46 even though they were not positive for other genes in either 
ST64B prophage. On sequencing of the 491 bp SB46 products it 
was found that they all had a sequence which had 9 SNPs 
compared with SB46 from ST64Box64- The sequence of SB46 
from ST64Bdtio4 has 3 SNPs (one located in one of the primers) 
compared with SB46 from ST64BrjT64 an d none of these is the 
same as the SNPs in the 9 SNP variant. When the SB46 product 
from a DT 1 04 isolate belonging to RGS was sequenced a mixture 
of the sequences from ST64B D xio4 an d the 9 SNP variant was 
found. BLASTing of the DT 1 04 genome sequence with the 9 SNP 
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SB46 sequence revealed a sequence with 100% identity to the 9 
SNP variant located between coordinates 2832676 and 2833166. 
This lies within the location of Gifsy-1 on the DTI 04 genome [22]. 
BLASTing of the downloaded Gifsy-1 sequence from DTI 04 
against the ST64B D T64 sequence showed extensive sharing of 
sequence between the two prophages including SB46 (detailed 
later). It was concluded that the isolates in RGs 3, 4 and 9 
probably also had a Gifsy-1 sequence like that in DTI 04. 
Supporting evidence for this is presented in results for Gifsy-1 
variants. The three closely related isolates with the apparent 
variant ST64B DT104 (see above) in RG7 (ATCC 13311, 
05P41568972 DTI and 09P84391818 DT41) were positive for 
the SB46 sequence with 9 SNPs. They were also shown later to 
have other features which are consistent with them having the 
Gifsy-l DT io4 prophage. 

The SB46 sequences in isolates which had tested positive for 
ST64B DT64 or ST64B DT10 4 were compared. SB46 for selected 
samples with each type of prophage was amplified and sequenced. 
Isolates with the ST64B D tio4 prophage had the DT104 sequence 
with 2 SNPs (excluding one in the primer) relative to SB46 in 
ST64B DT6 4 (AY055382). Isolates with the ST64B DT64 prophage 
in RG2A (mainly DT141) and in all the Repeats Groups in Major 
Group 2 with the exception of those in RG14 (DT8, DT9 and 
DT64) had a sequence with one SNP relative to SB46 in 
ST64B DT64 (AY055382). This is the same sequence for SB46 in 
ST64B DT64 in the sequenced strains SL1344, DT2, D23580, 
ATCC 14028 and CVM23701. Isolates in RG14 had a mixed 
SB46 sequence consisting of the same 9 SNP variant seen in 
isolates presumed to have the Gifsy-1 sequence found in DTI 04 
and either the same SB46 sequence seen in ST64B DX64 
(AY055382) or the one SNP sequence seen in all of the other 
Repeats Groups in Major Group 2. As will be detailed later the 
Gifsy-1 sequence in RG14 isolates appears to be a hybrid between 
Gifsy-1 in DTI 04 and that in LT2. There were differences in 
distribution of ST64B SB46 sequence variants for the different 
phage types in RG14 (Text S2). 

The SB46 gene sequence amplified from the presumed 
ST64Bdtio4 m isolates belonging to RG6B had 14 SNPs relative 
to the same sequence in ST64B DT10 4. The evidence indicates that 
the 14 SNP variant is located in the ST64B DT104 prophage rather 
than in the Gifsy-1 prophage and that the SB46 gene has been 
acquired from another serovar (Text S3). 

Comparison of Gifsy-1 sequences in sequenced strains 

We made a comparison of the downloaded Gifsy- 1 sequences 
from sequenced S. Typhimurium strains LT2, SL1344, DT104, 
DT2, D23580 and ATCC 14028. The sizes and genome 
coordinates for the six Gifsy-1 sequences are shown (Table 1). 
As a result of sequence comparisons (Text S4 and S5) we have 
been able to identify apparendy unique sequences on Gifsy- 1dtio4 
not shared by ST64B DT64 or ST64B D tio4 or any of the other 
Gifsy-1 prophages, and on Gifsy-l SL i3 44 and Gifsy-1 DT2 not 
shared by other Gifsy-1 prophages. These sequences were used to 
design primer pairs (Table 2) for PCR tests which could be used to 
determine what type of Gifsy-1 was present in S. Typhimurium 
isolates. There was almost no significant sequence on LT2 not 
found on one of the Gifsy- 1 sequences or on Gifsy-2 . However we 
located a sequence between 32761 and 35070 on Gifsy-lLT2 
(NCJ310392) which is not shared by Gifsy-2 or Gifsy-1 DT104 but is 
shared by Gifsy- Isl 1344 an d Gifsy-lDT2- This was used to select 
primers for detection of possible Gifsy- 1 lt2 m isolates which were 
negative for Gifsy- Isl 1344, Gifsy-lDT2 an d Gifsy-lDTi04- s P ecmc 
sequences. 



Table 1. Sizes and genome coordinates of Gifsy-1 prophages 



for sequenced strains of S. Typh 


murium. 






Strain 


Accession No. 


Size (bp) 


Coordinates 


LT2 ATCC 43971 


NC_003197 


48491 


2728552-2777042 


DTI 04 NCTC 13348 


NC_022569 


49233 


2796744-2845976 


DT2 (Sanger) 


NC_022544 


50748 


2719156-2769903 


D23580 


NC_016854 


50767 


2752645-280341 1 


ATCC 14028s 


NC_016856 


50772 


2780105-2830876 


SL1344 


NC_0 16810 


51187 


2726269-2777455 



doi:1 0.1 371 /joumal.pone.0086203.t001 



The unassembled sequenced strain Salmonella Serovar 
CVM23701 4,[5],12:i:- is located in RG12A (Figures SI and S2) 
but BLASTing showed that it lacks the SL1344 and DT2 Gifsy-1 
specific sequences. To locate the Gifsy-1 sequences in the contigs 
for CVM23701 4,[5],12:i:- the Salmonella enterica sequences in 
GenBank were BLASTed with the Gifsy-lDT2 sequence. This 
identified a sequence on the CVM23701 contig NZ_A 
BAO01000067 between coordinates 5611 and 9500 which is 
substituted for the Gifsy- 1dt2 specific sequences. BLASTing of 
Salmonella sequences with this substitute sequence showed that 
the first 1274 bp sequence was shared only with S. Javiana 
GA_MM04042433 but the following 1464 to 3890 bp of the 
sequence showed over 96% alignment with Gifsy-1 and Gifsy-2 
sequences of Typhimurium LT2 (as well as LT2-related strains) 
and with Gifsy-2 but not Gifsy-1 sequences for other Typhimur- 
ium strains. This part of the sequence was also found to contain a 
tandem repeat. In both Gifsy-1 and Gifsy-2 of LT2 it appears as 
three 22 bp repeats followed by an extended 39 bp repeat, two 
more 22 bp repeats and a 12 bp partial repeat. In the CVM23701 
Gifsy-2 the tandem repeat has the same structure but in the 
CVM23701 Gifsy-1 it has the repeats structure 22x3, 39, 22x3, 
39, 22x2, 12x1. In the D23580 Gifsy-2 the tandem repeat is like 
LT2 but lacks the 39 bp repeat. We designed primer pairs to be 
used for detection of Gifsy-lcvM237oi-specific sequences in the first 
1940 bases of the CVM23701 sequence as well as another pair 
from the later part of the sequence to amplify DNA shared by all 
Gifsy-2 sequences and the Gifsy-1 in LT2 (Table 2). In addition we 
chose primers to target the tandem repeat. Finally, a visual 
comparison of the five significandy different Gifsy-1 sequences 
together with Gifsy-3 from Typhimurium ATCC 14028S is shown 
(Figure 2). 

Testing of isolates for Sequences from Gifsy-1 variants 
Distribution of Gifsy-lDTi04-specific sequences. The 

three Gifsy- l D Tio4- s P ec ific sequences corresponded to positions 
29105-29948 (a hypothetical protein in Gifsy-1 phages of several 
other Serovars), 38294-39171 (Replication protein O in several 
others Serovars) and 31778-34537 (within the arfTfigene complex 
Accession No. AB 104436.1). The results of tests for the three 
Gifsy- loTHM-specific sequences in representative isolates from the 
panel of 2 1 4 well-characterised isolates (Figures SI and S2) showed 
that isolates in RGs 3, 4, 8, 9A and 9B which are all in Major 
Group 1 had all three sequences and also lacked the GipA 
sequence. These RGs included DT179, DT101, DTI 2a, DT104 
and most of the DTI 97 isolates. It was concluded that all the 
isolates in these RGs probably have the Gifsy-lDTi04 prophage. It 
was noted that all the tested DTI 01 and DTI 2a isolates had the 
Gifsy- 1 D t 104 sequences even though each phage type consisted of 
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Table 2. Primer sequences for phage loci, VNTR loci and CRISPR 1 and 2 


in Salmonella 


serovar Typhimuri 


um. 




Primer 


Sequence (5 -3 ) 


Product 
size bp 


Target 


Source 


DTI 04 GiflA-F 


CACTCCAGACCTGCGCGACGATA 


414 


Gifsy-1 DT104 


NC_022569 


DTI 04 GiflA-R 


CTGGTATGTGTGTCCCTCGCCG 








DTI 04 Gifl B-F 


CGGCACCAGAACGGATCACGGCT 


583 


Gifsy-1 DT104 


NC_022569 


DTI 04 Gifl B-R 


GTAAAGCCAGACGTCAGGGAAACCG 








DTI 04 GiflC-F 


CCAGATCGTCAACTGTGTCATTCG 


636 


Gifsy-1 DT , 0 4 


NC_022569 


DTI 04 Gifl C-R 


TTGTATGGGACCTACGTTGCC 








DTI 04 Gifl D-F 


CACAGAGCAACTGAATGCAATG 


340 


Gifsy-1 DT104 


NC_022569 


DTI 04 Gifl D-R 


GAGACGGATCTGGAACATCA 








DTI 04 Gifl E-F 


CGGCGGAAGGCTGTAGAAACTA 


644 


Gifsy-1 DTI 04 


NC_022569 


DTI 04 Gifl E-R 


CGTCGTCGTATTGACTCTTGGC 








SL1344 Gifl A-F 


CACGCCACTCACCGTTAGGC 


413 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gifl A-R 


TCGCCTGTGTAATGACTACTGG 








SL1344 Gifl B-F 


CCTCTTTGAACCCGGCTCGC 


441 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gifl B-R 


GCAGGATGGCGCATTCCGTG 








SL1344 GiflC-F 


CGATTAGCCGTGGTAGCAGGG 


495 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gifl C-R 


CTGAAGACCTAACTGCTGCCG 








SL1344 GiflD.l-F 


CTTCCAAGAGCAGCCCCTGTT 


419 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gifl D.I -R 


CGATAGAGCTCAAGCAGGCGC 








SL1344 Gifl D.2-F 


CCGGTTACCAGCCTTCATCAAG 


500 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gif1D.2-R 


CAAGATCGACATCACCTGGAAGG 








SL1344 Gifl D.3-F 


CCTTCCAGGTGATGTCGATCTTG 


560 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gif1D.3-R 


CAAGCGTGTTGCCATCCTTCAC 








SL1344 Gifl E-F 


GCGCTTACCACCTGTGTGGAA 


446 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gifl E-R 


CATAATCCTTACCGCGGCCGG 








SL1344 Gifl F-F 


CCGCGGCCTGCTGGTTAAG 


431 


Gifsy-1 SL1 344 


NC_016810 


SL1344 Gifl F-R 


CGTACCAGTCCGGATCGCG 








DT2 Gifl A-F 


CCCAAGTTAGAGATGCAAAACGCC 


517 


Gifsy-1 DT2 


NC_022544 


DT2 Gifl A-R 


CGTCGATCGACCATAATTTGGGC 








DT2 Gifl B-F 


CACAACCCGGCCATGCTCAGG 


421 


Gifsy-1 DT2 


NC_022544 


DT2 Gifl B-R 


ATAAACCCCCTAAACCCCCCG 








DT2 Gifl C.I -F 


GCAGGAACCACTCCTTGATAACG 


566 


Gifsy-1 DT2 


NC_022544 


DT2 Gifl C.I -R 


ATCCAGTTCCTTGAGACGAAGCG 








DT2 Gif1C.2-F 


CGCCAAGAGGCAACTAAGTATCTG 


386 


Gifsy-1 DT2 


NC_022544 


DT2 Gif1C.2-R 


GTGAAAAATGCACATAATGACAAGCG 








DT2 Gifl 24bp Rpt-F 


CGGGGGGGTTAGGGGGTTTAT 


244 


Gifsy-1 QT2 


NC_022544 


DT2 Gifl 24bp Rpt-R 


CGGTCAGCGAATGGTCATCTTTAC 








LT2 Gifl-F 


GTGCTGTATCCAGTAGAAGCC 


760 


Gifsy-1 LT2 


NC_003197 


LT2 Gifl -R 


TCAGTCAGACACTACCATCGC 








CVM23701 Gifl A-F 


GGCAGGAAGGGCTAATTCACGAGA 


768 


Gifsy-1 CVM23701 


ABAO0 1000067 


CVM23701 GiflA-R 


CAATAGAGGCCGTGGCACTTGC 








CVM23701 Gifl B-F 


GCAAGTGCCACGGCCTCTATTG 


534 


Gifsy-1 CVM23701 


ABAO0 1000067 


CVM23701 Gifl B-R 


GCGCATCACGACGAACGACG 








CVM23701 GiflC-F 


CGTCGTTCGTCGTGATGCGC 


569 


Gifsy-1 CVM23701 


ABAO0 1000067 


CVM23701 Gifl C-R 


CAGAACTTCGGTAGCGGATTTGC 








SPTR4-F 


ACCTGGATAAATGGGCTTATTCCAAGC 


474 


Gifsy-1 CVM23701 


ABAO0 1000067 


SPTR4-R 


CATTCAGCCATTCGGGAACAGGAATAC 








Gifsy-3-F 


GAGCTCAGCAACGTGTCGAAAGC 


734 


Gifsy-3,4028 


NC_016856 


Gifsy-3-R 


TTGGTCTGCGGGGACATCGCC 








DTI 04 SB52-F 


GCGGCGATGAAGGAATTCGG 


605 


ST64B DT104 


NC_022569 
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Table 2. Cont. 







Primer 


Sequence (5 -3 ) 


Product 
size bp 


Target 


Source 


DT104 ^R57-R 

YJ 1 1 \J*-r JDJZ r\ 


/Av_V~/A\JMVwrt/A\J/A\_\_/AVw\_/AVjVw\J 








DT1 04 ^R53-F 

YJ \ 1 U*+ JDJJ \ 


rAr n r n AGGrGAATrAArATGrAA 

v_/AVJ \J/A\JVJV_\J/Ar\ 1 V_nnV_n 1 U\_nn 


538 


J l OH-DDT104 


Kir 077569 

OZZJO_7 


DT104 ^R53-R 

YJ 1 1 \J*-r JDJJ r\ 


CTCC &C,C ACTC ACCC A&CiATG 

v_ 1 VwVwMVJV_r\v_ 1 Vwr\V_V_^_r\r\U/A 1 VJ 








DT1 04 <IR7fi 1 -F 
yj i i ut jdzo. i r 


CCGCC A A ATTrTf^TrTnf^f^TTATr 

LUULUnnn 1 ILIUILI Uuu I Inlv. 


264 


CTfiAR 

J 1 D*IDdtiq4 


Mr nn^&Q 

1 N V- UZ Z JOS 


DT104 <IR7fi 1-R1 
YJ I I Ut JDZO. I n I 


CCT A A A CTCCTTCCCCTCTC A ArT 

1 r\r\r\U 1 VJU 1 1 1 1 LnnL 1 








DT1 04 <IR7fi 7-F 


A A A CTTC A C A C A AIT A C CCTC C 


388 


j I D*IDdt104 


1 N V- UZ Z JU5 


DT104 ^R?ft 1 -R7 
L/ 1 1 Ut JDZO. 1 nZ 


C C A C CCTC ATTf^TrTf^ A A rTTT 

LLnLUL 1 Un 1 1 U 1 V_ 1 UnnL 1 1 1 


583 






DT104 <;R7fi ^-F1 

L/ 1 1 Ut JDZO. J r 1 


CC Af^TA 1 1 1 1 C&CCCC&CT&G 


553 


J I 0^007104 


Mr nji^&Q 

1 N V- UZ Z JU5 


DT1 04 ^R?ft 7-R 
YJ 1 lut JDZO.Z n 


rTArTf^rf^rTf^AA AATArTf^r; 

\— 1 r\\. 1 Vj\^UVw 1 Unnnn 1 r\\- 1 








DT104 ^R7fi 3-F 

L/ 1 1 U*T JDZO. j r 


cjcc a g g atg AGrTG r a cj a a r 


444 


J 1 OH-DDT104 


Mr 077569 

UZZJO-7 


DT1f)4 ^R7fi 3-R 

YJ 1 I U*T JDZO. J n 


CCACC hTCCTC a r Tr Trr tgtg 

v_v_/A\JVJfA 1 v_v_ 1 VJrAv_ 1 V 1 1 V_ 1 U 1 VJ 








DT1 04 ^R33-F 

YJ 1 I U*T JDJJ r 


GTrAGGrArTrATTAAAGGTTGrt^ 

VJ 1 V_r\vjvj\_r\\_ 1 \-r\ 1 1 /A/A/AVJVJ I I \JV_u 


424 


J 1 OH-DDT104 


Mr 077569 

INV_ UZZJO-7 


nT1 04 ^R^-R 

L/ 1 lUt JDJJ n 


CCTT &TCC&TCCTT &TC ACTCTCC 

\J\J 1 1 r\ 1 Lun 1 U\_ 1 1 r\ 1 V-rtVJ 1 \_ 1 UL 








DT104 ^R^^R-F 
l/ i lut jdj/jo r 


(^nrTr(^(^T(^(^TATmTATAA(^ 

\J\JV_ | LUU | VJVJ 1 r\ 1 1 1 r\ 1 r\r\\J 


359 


J 1 □'+DDT104 


Mr 077^ftQ 

1 N V_ UZZ JO _7 


DT104 ^R^7^R-R 

L/ 1 lUt JDJ/JO n 


CTTTC AC AC ATCCCT ATTTCTTC 

\_ 1 1 1 *_r\VJr\Utt 1 Ouv. 1 r\ 1 1 1^1 1 L_ 








DT1 04 ^RAQ-F 

YJ 1 lUt JDt?^r 


CJCC A A ACC ACCCCTTCCCTTTC 

\J 1 LUnnnVJLnLLUL 1 1 UUU 1 1 1 U 


426 


J 1 0^007104 


Mr 077^ftQ 

1 N V_ UZZ JO J 


DT1 04 ^R4Q-R 
YJ 1 lut JDt7 n 


T A A AC A CCCCTCCC ATC A A CTCCC 

1 r\r\/AV_r\UVJV-VJ 1 UV-V-M 1 VJ/AttO 1 OVJ 








DT104 ^R^O^I-F 

L/ 1 1 Ut JDJUJ 1 1 


CC ATCCCCCCCACJCCC AT ATTT 


381 


J 1 OHDQT104 


Mr 077^fiQ 

1 N V_ UZZ JO J 


DT104 ^R5051-R 

YJ \ \ o*t JDJUJ 1 r\ 


CC A ACCC A AC ACC ACTGTATGT A AC 








STTR1-F 


CACCACTACAACCCTCACCACCAT 


770 


45 bp TR 


Mr 003197 

1 NV_ UU J 1 _/ / 


STTR1-R 


CCC C C A CCCTT A C C C C C C C ATGT A 








STTR7-F 


CCCCC ACC C CTTC JC ACT 


594 


^0/4^/^0 hn TR 
J3/tj/ jU Up In 


Mr nn^ 1Q7 

1 N V_ UUJ 1 J / 


STTR7-R 


TCTTCC ACCCC AAACCT ATCT A 

1 VJ 1 1 LLnULULnnnUU 1 n 1 L 1 n 








STTR8-F 


A ATTA ATTf^rrf^f^ ATf^f^Tf^ A 

r\r\ 1 1 /AM 1 1 ULLUUn 1 UVJ 1 \Jr\ 


841 


1 0R/1 1 hn TR 
i yjoi \ \ o up i r\ 


Mr nn^ 1Q7 

1 N V_ UUJ 1 J / 


STTR8-R 


ACCCATTCCTCCCCT ACAT 

MULUn 1 1 UL 1 UULL 1 nun 1 








STTR1 1-F 


ATCCAACCCCTCCTTACCTC 


71 6 


1 hn TR 
i jj up i n 


Mr nn^ 1Q7 

1 N V_ UUJ 1 J / 


STTR11-R 


ACGTAGCCCCGTATCTGATG 








STTR12-F 


GGATGGTGGTGTTATTGGCCGGT 


762 


184 bp TR 


NC_003197 


STTR12-R 


CTGAAGGAGCACGCCTGGAAAGTG 








CDC ST3-F 


GTTCTTCTGCAACGCAGGCA 


193 


12 bp TR 


NC_003197 


CDC ST3-R 


GATGGCATGACGCTGCAACG 








CRSP1-F 


CGAAGGCGGAAAAAACGTCCTG 


1735 


CRISPR1 


NC_003197 


CRSP1-R 


GACGTATTCCGGTAGATCTGGATG 








CRSP1 74 bpSp-F 


CGCGGGGAACACAATTAAAGCCGA 


1 273 


CRISPR1 after 74 bp 
spacer 


NC 003197 


CRSP1 74 bpSp-R 


TCGGCTTTAATTGTGTTCCCCGCG 


486 


CRISPR1 before 74 bp 
spacer 




CRSP2-F 


GTASC WGCC ATTACTGGTACAC AG 


2197 


CRISPR2 


NC_003197 


CRSP2-R 


CATAGCGATGCACGGATCACGC 








CRSP2-Sp13-F 


ATTTCGCCTTCGGCACTGACGTCAC 


1396 


CRISPR2 after spacer 13 


NC_003197 


CRSP2-Sp13-R 


GTGACGTCAGTGCCGAAGGCGAAAT 


827 


CRISPR2 before spacer 
13 


NC_003197 


CRSP2-Sp24-F 


CTAG G AGG CGTA ATG AATACTACG 


400 


CRISPR2 after spacer 24 ABAO01 000021 


CRSP2-Sp24-R 


CGTAGTATTC ATTACG CCTCCTAG 


1434 


CRISPR2 before spacer 
24 


ABAO0 1000021 



doi:1 0.1 371 /journal.pone.0086203.t002 



isolates with distinctly different genotypes (Figures SI and S2). 
Isolates in RG14 in Major Group 2 consisting of DT8, DT9 
and DT64 isolates had the GipA sequence and two of the Gifsy- 
1dtio4 sequences but lacked the sequence corresponding to the 



29105-29948 insert in Gifsy-lDTio-t- These isolates appear to have 
a hybrid form of Gifsy- 1 in which the first part of the genome to 
beyond the insertion point for the insert is like that seen in 
LT2, DT2 or SL1344 and the remainder is like that seen in 
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Salmonella enterlca subsoenierica 14028s Gifsy-3 



Figure 2. Comparison of Gifsy-1 sequences for 5. Typhimurium strains LT2, DTI 04, SL1 344, DT2 and Salmonella Serovar CVM23701 
and Gifsy-3 sequence for 5. Typhimurium ATCC 14028s by MAUVE (See also Text S5 for further explanation). 

doi:1 0.1 371 /journal.pone.0086203.g002 



Gifsy-1 dti 04- We have designated this prophage Gifsy- 1dt64- 
Sequencing of a selection of PCR products for all three Gifsy- 
lDTi04-sp ec ific sequences from each of the RGs showed no 
difference from the sequences in DTI 04 except for isolate 
06P23331584 in RG3 (Figure SI) (Text S6). 

Distribution of Gifsy-lsLi344-specific sequences. The 
results of tests for the six Gifsy-1 sLi344-specific sequences in 
representative isolates from the panel of 214 well-characterised 
isolates (Figures SI and S2) showed that isolates in RG10A and 
RG10B (mostly DT44, DT29 and DT35 isolates, along with some 
DT120, DT186, U302 isolates and including strain SL1344) in 
Major Group 2 had all six sequences as well as the GipA sequence 
and therefore were presumed to have the Gifsy-lsLi344 prophage. 
The DT126 isolates in RG11A and in RG11B (also in Major 
Group 2) had all six Gifsy-1 SL1344 sequences as well as the GipA 
sequence but subsequent testing showed that they lacked the Gifsy- 
1lt2 sequence also found in Gifsy- Isli 344 and Gifsy-1 dt2 but not 
in Gifsy-1 DT1{)4 . This result was supported by the same result in a 
further 19 DTI 26 isolates. As this indicates a significantly different 
sequence from Gifsy-lsu344 we have designated this prophage 
Gifsy-l D Ti26- I n Major Group 1 RG2 isolates which were mostly 
phage typed as DT 1 4 1 and variants of DT 1 4 1 but not any of the 
DT4 isolates (including LT2) had all six of the Gifsy- Isli 344 
sequences as well as GipA and the Gifsy-1 LT2 sequence also found 
in Gifsy-lsu344 and Gifsy-lnTz but not in Gifsy-lDTi04- They 
were presumed to have the Gifsy-lsLi344 prophage. All the DT141 
and DTI 41 variants which were tested for CRISPR sequences had 
the RG2A-type CRISPR 2 sequences. 

Distribution of Gifsy-lDT2 _s P ec if> c sequences. The re- 
sults of tests for the three Gifsy-1 DT2 _ specific sequences in 
representative isolates from the panel of 214 well-characterised 
isolates (Figures SI and S2) showed that isolates in RG1 (mainly 
DTI 70 and DTI 2 but with some DTI 86 and DTI 93 isolates) in 
Major Group 1 had all three sequences as well as GipA and were 



presumed to have the Gifsy-1 dt2 prophage. Selected isolates in all 
the divisions of RG12 and RG13 in Major Group 2 as well as 
06P75527334 in RG10A (Figure S2) also had all three sequences 
as well as GipA and were presumed to have the Gifsy-1 dt2 
prophage. RG12 and RG13 were populated with mostly DT135, 
DT135a* and DT170 isolates with smaller numbers of U307, 
DT102, DT6 var 1, U290 and DT3 isolates. Significantly, isolate 
06P75527334 was typed as DTI 35a which was different from the 
other phage types in RG10A. Its location there was determined 
solely by its CRISPR 2 profile since its CRISPR 1 profile placed it 
in RG12D (Figure SI). Isolate 06P75527334 is one of a number of 
hybrid isolates in the panel in which the CRISPR 1 and 2 profiles 
belong to different Repeats Groups. The sequenced strains ATCC 
14028 and D23580 have been placed in RG12A and RG12B 
respectively and have the Gifsy-1 dt2 prophage but DT2 itself has 
not been assigned to a Repeats Group because its CRISPR profile 
does not conform to any of the recognised Repeats Groups. Isolate 
08P8821 1 137 in RG12A was found to have all three Gifsy- 1 DT2 - 
specific sequences. This strain had previously been found to have 
all six Gifsy-lsLi344-specific sequences. It was subsequently shown 
that these sequences are located on the Gifsy-3 prophage present 
in this strain (see below). One isolate assigned to RG2B had all 
three Gifsy- 1 DT2" s pecific sequences and another also in RG2B had 
just one Gifsy- l D T2" s pecific sequence but as will be shown later 
most RG2B isolates have the Gifsy- 1lt2 variant. 

Identification of Gifsy-lsequences specific to CVM23701 
in selected isolates. Using the primer pairs derived from 
CVM23701 to test a limited range of isolates we showed that the 
CVM23701-specific sequences and the larger tandem repeat 
sequence found in this strain were present in only one local isolate 



*DT135a is a variant of DT135 recognised by Australian Salmonella phage typing 
laboratories but not given official status by the Health Protection Agency located 
in Colindale, London. 
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(07P69570074) from 2007 which had a RDNC phage typing 
result. This isolate significantly lacks SL1344 and DT2 Gifsy-1 
specific sequences and is diphasic not monophasic. This isolate has 
the same MLVA and CRISPR profiles as CVM23701 and 
appears to have the same Gifsy-1 prophage as CVM23701. We 
have therefore designated this prophage Gifsy-1 cvM2370i- We 
obtained DNA from a strain of DTI 9 la isolated in the UK in 
2008 and associated with an outbreak of human salmonellosis 
traced to handling of feeder mice imported from the USA [24,25]. 
Like CVM23701 this strain is monophasic. It also lacks the 
SL1344 and DT2 Gifsy-1 specific sequences but has the two 
CVM23701 Gifsy-1 -specific sequences. However it lacks the larger 
tandem repeat sequence found in CVM23701 and local isolate 
07P69570074. All isolates tested for the tandem repeat located on 
Gifsy-2 (and Gifsy-1 in LT2) had the same size allele and not the 
shorter allele found in D23580. 

Identification of isolates likely to have the Gifsy-1 LT2 
variant. Since there were no unique sequences found in Gifsy- 
1 LT2 we applied a PGR test for amplification of a Gifsy-1 lt2 
sequence found also in Gifsy-l S Li344 an d Gifsy-lDT2 but not in 
Gifsy-2 or Gifsy-1 dtkm as a proxy test for the possible presence of 
Gifsy-1 L T2- Isolates negative for Gifsy-lDTi04, Gifsy-lsLi344 an d 
Gifsy-1 D T2 -specific sequences but positive for GipA and STTR6 
were tested. Isolates in RG2 (DT4, DTI 93) which had not tested 
positive to one of the other Gifsy-1 types as well as additional DT4 
isolates were all positive for the Gifsy-1 LX2 sequence. They most 
likely have Gifsy-1 LT2 because LT2 is DT4 and belongs to this 
Repeats Group. All tested isolates with the 171-x-x-0-489 MLVA 
in RG6B were positive as were all DT 1 3 isolates, a related DT4 1 
isolate and an anomalous DT197 isolate in RG1 1A as well as a 
DTI 60 isolate and a closely related untypable isolate, the only 
members of RG15. Sequencing of the Gifsy-1 lt2 sequence from a 
DT4 and a DT193 (RG2) isolate showed 100% identity with the 
sequence in LT2. However, sequencing of the same product from 
four isolates with the 1 7 l-x-x-0-489 MLVA belonging to RG6B as 
well as the untypable strain 08P662 17254 (in RG1 in Figure SI 
and RG4A in Figure S2) and 07P25565432 typed as DT197 (but 
with a unique STTR1 allele and CRISPR 1 and 2 sequences most 
closely related to those of DTI 26 isolates in RG11A) showed 
alleles with a single SNP at base 643. 

As a validation test we tested a range of isolates positive for 
Gifsy-lrjTi04, Gifsy-lDT64j Gifsy-lsLi344 or Gifsy-lDT2 sequences 
for the presence of the Gifsy-1 lt2 sequence. All of the isolates 
positive for Gifsy-lDTi04 or Gifsy-lDT64 were negative for the 
Gifsy-1 lt2 sequence. All of the isolates from various Repeats 
Groups which were positive for Gifsy-1 sli 344 or Gifsy-1 dt2 were 
positive for the Gifsy-1 L T2 sequence except for DTI 26 and DTI 26 
var 4 isolates in RG11A and B, which had previously tested 
positive for all of the Gifsy-1 sli 344 sequences and yet tested 
negative for the Gifsy-1 lt2 sequence. As previously indicated we 
have designated this prophage Gifsy-1 D T126- 

Distribution of Gifsy-3 

14028 sequence 

Isolate 08P88211137 in RG12A, phage typed as DT135, was 
positive for all six Gifsy-1 S l 1344 sequences even though other 
DTI 35 isolates were negative for all six. It was found that this 
isolate was also positive for Gifsy- 1 DT2~ s P ec ific sequences as were 
other DTI 35 isolates. This isolate is one of a small group of 
isolates placed in RG12A of Major Group 2 in Figures SI and S2 
because the CRISPR 2 profile was the same as for strains 
CVM23701 4,[5],12;i;- and ATCC 14028s although the CRISPR 
1 profiles were the same as other DTI 35 and DTI 02 isolates in 
RG12B in Figure SI. BLASTing of the ATCC 14028s sequence 
with the Gifsy-1 SL i 344 sequence showed that there was an 



alignment at 99% identity but in reverse order for the Gifsy- 
1sli344 sequence 36491 to 46388 (Figure 2) with ATCC 14028s 
sequence (NC_0 16856) at coordinates 1289242 to 1299140. This 
included the whole region where the six Gifsy-lsu344 sequences 
are located. This was known to be the approximate location of the 
Gifsy-3 prophage [21]. Following this discovery it was noted that 
the MLVA profile of isolate 08P8821 1 137 (DT135), 162-300-318- 
357-523, was very similar to that for ATCC 14028s, 162-306-318- 
357-523. We therefore explored our S. Typhimurium isolate 
database from 2006 onwards for MLVA profiles the same as those 
for 08P8821 1 137 or ATCC 14028. We found seven isolates with 
the same MLVA as 08P882 11137 including two from kangaroo 
meat and one with the same MLVA as ATCC 14028s. All were 
typed as DTI 35 except one of the isolates with the 08P8821 1 137 
MLVA which was RDNC and the two kangaroo isolates one 
typed as DTI 2 and the other as RDNC. On testing for Gifsy-1 
sequences all were positive for both Gifsy-1 sli 344 an d Gifsy- 1dt2 
sequences. By alignment of the downloaded 50884 bp Gifsy-3i4028 
sequence (ATCC 14028s genome coordinates 1284562 and 
1335445) with the various Gifsy-1 sequences and the downloaded 
Gifsy-2 sequence we located a sequence at 46191 to 47760 on 
Gifsy-3i4o 2 8 (Figure 2) with genome coordinates 1330752 and 
1332321 which appeared to be specific to Gifsy-3i4( )2 8. We chose 
primers for a 734 bp internal sequence and tested all of the isolates 
with MLVA same or similar to ATCC 14028s and a represen- 
tative sample of 33 isolates positive only for the Gifsy-lsu344 - 
specific sequences. Only the seven isolates previously positive for 
both Gifsy- Isli 344 and Gifsy- 1dt2 sequences were positive for the 
Gifsy-3i4o 28 sequence. PCR product from two of the positive 
isolates was sequenced and shown to have 100% identity with the 
Gifsy-3 140 28 sequence. 

Correlation of Prophage Profiles with Repeats Groups 
and Phage Types 

On the basis of the two ST64B prophages including two 
variants of ST64B D tio4, seven Gifsy-1 prophages, Gifsy-3 
prophage and a variant of SB26 gene from ST64B DX104 we were 
able to distribute almost all S. Typhimurium isolates into sixteen 
prophage profiles (Table 3). Some profiles were exclusive to one 
Repeats Group while others were common to two or several 
Repeats Groups. RG11A contained two profiles and RG12A 
contained three. There was a high degree of correlation between 
prophage profiles and phage types with many phage types falling 
completely or almost completely within a single prophage profile. 
Some phage types, notably DTI 20, DTI 70, DTI 86, DTI 93, DT6 
var 1, DTI and U302, were represented in two or more prophage 
profiles usually consistent with the Repeats Group to which they 
had been assigned. Thus DTI 70 isolates fell almost entirely into 
two distinct Repeats Groups (1 and 13 A) which correlated 
completely with prophage profiles 6 and 8 respectively. Both 
profiles had Gifsy- 1dt2 but profile 6 isolates had ST64B D xio4 ar >d 
profile 8 isolates had ST64B DT64 . 

Discussion 

Previous studies have reported significant sequence differences 
for ST64B and Gifsy-1 prophages from different strains of S. 
Typhimurium but they have by and large not investigated these 
differences further and they have not mapped the distribution of 
these variants among isolates with a diverse range of genotypes 
and phage types. By comparison of the available sequenced 
prophages we have shown that there are two significant variants of 
ST64B: the prototype from DT64, which is almost the same in 
strains SL1344, DT2, D23580, ATCC 14028s and CVM23701 
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4,[5],12:i:-; and the one in DT104. We have shown likewise that 
there are five significant variants of Gifsy-1 amongst sequenced 
strains, one in LT2, one in DT104, one in SL1344, one in 
CVM23701 4,[5],12:i:- and one shared by DT2, D23580 and 
ATCC 14028s. By comparison of the variants of each prophage 
we have identified sequences which appear unique to a particular 
prophage variant or, in the case of Gifsy-1 lt2> shared by a number 
of variants. Our investigations have led to the rediscovery that 
significant gene sequences are shared between ST64B and Gifsy- 1 
prophage from DTI 04 and made the further rediscovery that a 
large segment of Gifsy-1 from SL1344 is shared by Gifsy-3 from 
ATCC 14028s. We have used the techniques of PCR and 
sequencing of PCR products to test for the presence of these 
'unique' sequences in S. Typhimurium isolates, from a large 
collection of local isolates, chosen to represent most of the phage 
types encountered locally and most of the genotypes we have 
identified by VNTR and CRISPR analysis over the period 2000 to 
2012. As a result of extensive testing we have identified two new 
variants of Gifsy-1, one in DT126 isolates and one in DT8/9/64 
isolates, as well as two variants of ST64B D xio4- We have now 
found that the sequenced S. Typhimurium strains UK-1 
(NCJ316863) [26] and 08-1736 (NC_021820) contain a Gifsy-1 
prophage which is very closely related to the one in CVM23701. It 
has the same sequence which we have shown to be unique to 
CVM23701 and a local strain as well as the extended 22 bp 
VNTR also only found in those two strains. 

To identify a particular prophage variant in an isolate we have 
made the assumption that an isolate which has all of the sequences 
thought to be unique to a particular variant as well as all of the 
tested sequences shared by that variant with other variants is likely 
to have a prophage which is the same as or closely related to that 
variant. Support for this contention comes from the S. Typhimur- 
ium whole genome sequences deposited as Whole Genome 
Shotgun projects at DDBJ/EMBL/GenBank under the accessions 
AMDX, AMDY, AMEA, AMEB, AMEC, AMED, AMEE, 
AMEF, AMEG, and AMEH. All the isolates are strains of 
DTI 35a and belong to the RG13A genotype. The Gifsy-1 
sequences on most of these strains have been located and found 
to be almost identical to the Gifsy- 1dt2 variant sequence with the 
notable exception of a 24 bp VNTR located at 2172456 to 
2172555 on the assembled DT2 genome (Sanger) (Text S7). The 
sequence for S. Typhimurium strain TN061786 under accession 
AERV0 1000000.1 has a Gifsy-1 sequence almost identical to the 
Gifsy- 1dt2 variant sequence as well as the same 24 bp VNTR. 
TN061786 is a hybrid strain which like local isolate 08P662 17254 
belongs to RG1 in Figure SI and RG4 in Figure S2. Further 
support from another angle comes from the Southern hybridisa- 
tion analyses of Tucker and Heuzenroeder [17] (Text S8). 

We have been able to map the distribution of Gifsy-3 and 
variants of ST64B and Gifsy-1 prophages among local isolates. We 
have found a very high degree of correlation between the VNTR/ 
CRISPR genotype which we have termed Repeats Group and the 
prophage profile based on ST64B and Gifsy-1 variants (Table 3). 
We believe that this correlation lends support to the concept of 
VNTR/CRISPR genotypes or Repeats Groups in S. Typhimur- 
ium. There is a very similar system of genotyping applied to 
isolates of Mycobacterium [27,28,29]. Some Repeats Groups such as 
RGs 1 , 6B, 8 and 1 4 have all or nearly all members with the same 
unique prophage profile. Others may share their prophage profile 
with one or more Repeats Groups while some Repeats Groups are 
divided into two or more different profiles. Further differentiation 
can be obtained by including Gifsy-3 for some members of 
RG12A and the SB26 DT i 0 4 variant found in DTs 197, 43 and 6 in 
RG9A. Results of phage gene profiling [19] show that more 



differentiation can be obtained by including markers from other 
phage groups such as P22- and P2-like phages. 

Some phage types such as DTs 4, 8, 9, 12a, 22, 44, 64, 101, 179 
and 197 correlate totally or almost so with a single prophage 
profile even for DTs 12a and 101 which are divided into two 
different Repeats Groups. Others such as DTs 170 and 6 var 1 
correlate with two distinct prophage profiles. Thus isolates with 
one prophage profile have a genotype and a Repeats Group which 
are different from isolates with the other prophage profile. Still 
other phage types such as DTs 120, 186 and 193 have quite a 
diverse range of genotypes and prophage profiles. This is most 
likely because these phage types show very few reactions or none 
(DTI 93) with the Anderson typing phage panel. 

We have tried to assess how much influence the different 
variants of ST64B and Gifsy-1 may have on phage typing. In the 
case of DTs 4 and 1 4 1 there is evidence that the correlation of 
DT4 with the Gifsy- 1 L T2 variant and of DT141 with the Gifsy- 
l SL i344 variant may be related to the difference in phage typing 
pattern. It may be supposed that since most of the Anderson 
typing phages are P22 phages [30] it would be likely that phage 
type would be more greatly influenced by the type of P22 phage in 
an isolate. It is known that conversion of DT9 to DT64 and of 
DTI 35 to DTI 6 is caused by lysogenisation with ST64T, one of 
the P22 phages [3]. However DT4 and DT141 lack P22 phages by 
phage gene profiling [19] so it seems possible that the different 
Gifsy-1 variants may be contributing to the observed phage type 
difference. Arguing against this is the example of isolates with the 
phage type 6 var 1. Some have the combination of ST64BrjT64 
with Gifsy- 1dt2 while others have the combination of the 
SB26dtio4 variant with Gifsy-lDTi04 but no ST64B prophage. 
There appears to be no P22 phage by phage gene profiling so the 
phage type appears to be unaffected by the difference in prophage 
profiles. Perhaps some Gifsy-1 variants influence phage typing 
patterns but others do not. The example of DTI 70 where isolates 
of two distinct genotypes with prophage profiles which have 
different ST64B variants but the same Gifsy- 1 and apparently the 
same P22 prophage [19] produce the same phage typing pattern 
appears to indicate that the difference in ST64B phage sequence 
does not affect the phage type. Some DT141 isolates have 
ST64BrjT64 while others lack the prophage and yet the phage 
typing result is not affected. However, Tucker and Heuzenroeder 
[17] observed changes in phage type following introduction of the 
imm-C-\\kt genes from ST64B DT64 into isolates. 

We have divided all S. Typhimurium isolates into two groups 
termed Major Groups 1 and 2. The division is based mainly on the 
different alleles of STTR3, STTR7 and STTR9 (Text S9). We 
have noted that isolates in Major Group 2 practically always have 
the SB26 gene from ST64B DTB4 while isolates in Major Group 1 
always lack the gene with the exception of a minority of members 
of RG2 A, mostly DT 1 4 1 . We also note that the phage profiling 
results indicate that ST64B prophages seem only to occur when a 
Gifsy-1 prophage is present in an isolate. Thus the marked 
difference in occurrence of ST64B DX64 between the two Major 
Groups is possibly related to the type of Gifsy-1 variant or 
combination of Gifsy- 1 variant and ST64B variant present in an 
isolate. All isolates in Major Group 2 have ST64B D T64 combined 
with one of six Gifsy-1 variants. In contrast isolates in Major 
Group 1 may have both ST64B and Gifsy-1 prophages, Gifsy-1 
but no ST64B or lack both ST64B and Gifsy-1 prophages. Four of 
the Repeats Groups in Major Group 1 have ST64B D xio4 or one of 
its variants in combination with one of three Gifsy- 1 variants and 
RG2A has some DT141 isolates with ST64BrjT64 combined with 
Gifsy- l SL i 344. The latter combination is the only one which also 
occurs in Major Group 2 isolates. Gifsy- 1dt2 occurs in both Major 
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Groups but in Major Group 1 it is combined with ST64B D tio4 
while in Major Group 2 it is combined with ST64B DX64 . Apparently 
Gifsy-lLT2 occurs in both Major Groups but only in RG15 isolates 
(DT 1 60 and related Untypable isolate) and DT 1 3 isolates belonging 
to RG1 1 A is it accompanied by ST64BrjT64- Since it is known that 
the phage type of a strain can be altered by the introduction of a new 
phage lysogen it is not surprising that there is a degree of correlation 
between phage type and prophage profile. However it is difficult to 
explain the even higher level of correlation between the repeats 
profile and prophage profile even though the STTR6 VNTR is 
located on Gifsy-1, and CRISPR sequences are known to carry 
spacer sequences which can be derived from prophages and are 
known to confer immunity from infection with phages which carry 
the same sequences as in the spacers [31]. The composition of the 
Repeats Groups depends as much on the combinations of VNTR 
alleles as it does on the CRISPR profiles. 

It is possible that superinfection exclusion systems controlled by 
residing prophages is driving these divisions. Gifsy-lDTi04 an d 
probably Gifsy-1 DT64 , from the results of Tucker and Heuzenroe- 
der [17], have the same immC region (SB38, 39 and 40 in the 
ST64B DT64 genome) as ST64B D T64 an d ST64B D tio4 but the 
other Gifsy-1 variants do not. Consequendy it may be expected 
that any isolate with Gifsy-lDTi04 or Gifsy-lDT64 ma y be immune 
to infection by an ST64B phage. There is partial evidence for this 
in the table of prophage profiles. Thus four out of six Repeats 
Groups in Major Group 1 with Gifsy-l D Ti04 d° n °t have any 
ST64B prophage. However, in the case of isolates in RG8, mainly 
DT104, dual lysogeny with ST64B D tio4 and Gifsy-lDTi04 is 
permitted. It should be noted that co-existence of two prophages 
with the same immunity module occurs in LT2 with Gifsy-1 and 
Gifsy-2 [15]. It is presumed that this is the result of an 
intrachromosomal recombination/conversion event [13]. This 
may be the case also for ATCC 13311 and related members of 
RG7. However the ST64Brjxio4 variant in these isolates may not 
have the ST64B immC region because it has the two ST64Bqtio4" 
specific sequences before the immC region but not the ones after it. 
All the isolates in Major Group 2 have ST64B DT64 even in RG14 
isolates which have Gifsy-1 DT64 which possibly has the same 
immunity module as ST64Box64- RG14 isolates may therefore 
represent another example of recombination producing lysogeny 
with two prophages with the same immunity module. Even 
allowing for exceptional recombination events it is still possible that 
superinfection exclusion systems are responsible for maintaining 
clones and hence for the strong correlations we see between 
genotype, phage type and prophage profiles. For example, it is 
curious that members of RG1 in Major Group 1 with Gifsy-lryr2 
have ST64B DX1()4 while all the isolates in Major Group 2 with 
Gifsy-1 dt2 have ST64B DT64 . Then again, members of RG2B in 
Major Group 1 with Gifsy-1 SL1344 ma y or may not have 
ST64BrjT64 while all members of RG 1 0 in Major Group 2 with 
Gifsy-1 SL1344 have ST64BrjT64- It may be that certain combinations 
of prophage are difficult to disrupt and so clones arise because they 
are protected from the genetic changes resulting from invading 
phages. These clones survive if they have selective advantage but are 
gradually replaced as new strains with a more advantageous mix of 
prophages and other mobile elements arise. This is consistent with 
the conclusion of Figueroa-Bossi et al 2001 [12] that lysogenic 
conversion is a major mechanism driving the evolution of Salmonella 
bacteria. There are also P22 and P2-like prophages which if present 
are likely to contribute to the maintenance of clones. 

We have identified an emerging S. Typhimurium genotype with 
phage types DT193, DT120, DT195, U302, RDNC and 
Untypable with a 1 7 1 -x-x-0-489 MLVA profile and multiresistant 
antibiotic profiles. Possibly the same genotype has been reported 



in the United Kingdom [29] where DTI 20 isolates had similar or 
same MLVA and antibiotic resistance profiles as seen here. A 
study of Salmonella enterica Serovar 4,[5],12:i:- monophasic variants 
of S. Typhimurium isolated in Germany from pigs, pork and 
humans in 2006 and 2007 [32] showed that the most common 
phage types were DTI 93 and DTI 20, most MLVA profiles were 
1 7 1 -x-x-0-489 and most were multiresistant with the most 
common resistance profile the same as reported in our isolates 
and in the UK. It would be interesting to compare the prophage 
profiles for the isolates from both countries to see if they are the 
same as for the emerging genotypes we are seeing in Australia. It 
should be noted that strain CVM23701 4,[5],12:i:- included in our 
table of isolates is a different genotype from nearly all of the 
isolates in the German study and it has a different prophage profile 
from our isolates with the 1 7 1 -x-x-0-489 MLVA (see [33]). 

We have made the interesting discovery that isolates in RG9A 
with phage type DT197, DT6 and DT43 have a sequence variant 
of the ST64B D tio4 gene SB26dtio4 but have no other tested 
ST64B D tio4 sequences. Isolates with phage type DT12a (RG9B) 
lack the SB26dtio4 variant. They are placed in RG9B because the 
MLVA profiles are closely related to those of DTs 197, 6 and 43, 
the CRISPR 1 profiles are identical to these other phage types and 
the CRISPR 2 profiles are related, but DTI 2a isolates have 
significandy more spacers than DTs 197, 6 and 43 isolates. It 
seems probable that SB26dtio4 is located on another prophage. It 
also seems likely that DTs 197, 6 and 43 have arisen from DT12a 
(RG9B) by acquisition of further prophage or prophages and a 
simultaneous deletion of CRISPR 2 spacers. We note from our 
phage gene profiling results [19] that DTI 2a isolates have markers 
for a P22 prophage not present in DTs 197, 6 or 43 isolates but 
DTs 197, 6 and 43 isolates have the SopE marker from SopEfl), a 
P2-like prophage, not present in DTI 2a. However testing (not 
shown) indicates that the SopEO insertion point in DTI 97 isolates 
is vacant and the SopE product from DTI 97 isolates has 5 SNPs 
relative to the SopE sequence in SopEfl). Therefore the SopE gene 
variant in DTI 97 must have another location perhaps in another 
phage. Sequencing of DT197 and DT12a strains would be 
desirable to elucidate their relationship. 

The detection of seven isolates among thousands with the rare 
Gifsy-3 prophage by the similarity of their MLVA profiles to that 
of ATCC 14028 is strong evidence that MLVA profiles have a 
significant probability of predicting phylogenetic relationships 
between strains. It also indicates that clones of related strains are 
relatively constant over time since ATCC 14028 was isolated in 
1960 [21]. Clonal constancy is also demonstrated by the 
unexpected discovery of relatedness between ATCC 13311 and 
two local isolates from 2005 and 2009 which have almost the same 
MLVA as ATCC 13311 and the same unique ST64B DT104 variant 
as well as Gifsy-1 DT104 . ATCC 13311 was isolated in 1911 
probably in the Netherlands [34]. 05P4 1568972 DTI was isolated 
in 2005 from an adult male who is the neighbour in a suburb of 
Brisbane of the adult male from whom 09P84391818 DT41 was 
isolated in 2009. The latter patient is thought to have acquired his 
Salmonella infection in Vietnam. 

The two members of RG15 are: MDU 2009 16 which was 
isolated in Tasmania and typed as DTI 60 which is the same as the 
phage type associated with infection in birds and humans in New 
Zealand since 2001 [35]; and 05P08494937 which was isolated in 
2005 from a Queensland patient and was phage Untypable. The 
MLVA profiles for the two isolates are closely related. Both have 
unique STTR3, STTR7 and STTR8 alleles and an STTR12 
allele rarely seen in isolates other than DTI 26 and DTI 3 isolates. 
Both CRISPR 1 and 2 profiles are unique. The CRISPR 1 profiles 
are the same for the two isolates but the CRISPR 2 profiles are 
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significantly different although clearly related. Both isolates have 
ST64B D T64 an d probably Gifsy-lLT2 as f ar as ° ur tests can reveal. 
This prophage combination is only shared with DTI 3 isolates in 
RG11A. Clearly there is potential for this type of geno typing to 
discover more diversity among S. Typhimurium isolates from 
other animal and environmental niches. 

Conclusions 

As a result of a survey of characterised S. Typhimurium isolates 
seeking specific sequences of prophage variants we have discov- 
ered at least two additional variants of Gifsy-1, one in isolates of 
DTs 8, 9 and 64 which appears to be a hybrid of non-DT104 
Gifsy-1 and Gifsy-1 dtk)4, another which appears to be specific to 
DTI 26 isolates and recognised a third in Salmonella CVM23701 
(GenBank ABAO00000000. 1) and a local isolate. We have also 
made the rediscovery that Gifsy-1 in SL1344 shares specific 
sequences with Gifsy-3 in ATCC 14028. A search for isolates 
which have an MLVA profile the same as or similar to that of 
ATCC 14028 has allowed us to identify, among several thousands, 
a few rare isolates which harbour this prophage. On the basis of 
two ST64B prophages including two variants of ST64B D tio4> 
seven Gifsy- 1 prophages, Gifsy-3 prophage and a variant of the 
substituted SB26 gene from ST64B D tio4 we were able to divide 
almost all S. Typhimurium isolates into sixteen profile groups. 
Each of these groups is specific, or largely so, to a range of phage 
types. In addition we find that each prophage profile is specific to 
one or a number of genotypes (Repeats Groups) as determined by 
the combination of VNTR and CRISPR typing. 

Materials and Methods 

Bacterial isolates and DNA preparation 

The isolates of S. Typhimurium used in this study were nearly all 
obtained from specimens collected in Queensland and northern 
New South Wales. They were mostiy derived from cases of human 
infection. The majority of them were locally-acquired but significant 
numbers could have been acquired overseas. They were recruited 
mainly from the approximately 2,500 isolates routinely typed by 
MLVA at the Public Health Microbiology Laboratory in the period 
2006 to 2009 and from a similar number of retrospectively MLVA- 
typed isolates from the 1999-2005 period. A panel of 214 isolates 
which had been well characterised and assembled into related 
groups termed Repeats Groups (RGs) according to their VNTR and 
CRISPR profiles was tested most intensively (Figures SI and S2). 
When correlations between phage gene test results and RG or other 
groups of isolates were observed, further isolates with same or similar 
MLVA profiles or phage type were selected to test the correlations 
more rigorously. In addition, a culture of S. Typhimurium ATCC 
13311 was obtained from the American Type Culture Collection. 
DNA was prepared by taking one colony from an overnight culture 
into 400 ju.1 of TE buffer and boiling for 8-10 minutes. 

Comparison of ST64B, Gifsy-1 and Gifsy-3 Sequences 

The sequence for ST64B DX64 Accession No. AY055382 and the 
gene map were obtained from GenBank (http:/ /www.ncbi.nlm. 
nih.gov/genbank/). The whole genome sequences for S. Typhi- 
murium DT2 and DTI 04 were downloaded from the Sanger 
website (http://www.sanger.ac.uk/Projects/Salmonella). The 
Align program in NCBI website was used to align the ST64B DX64 
sequence with each of the DT2 and DTI 04 sequences to locate 
the boundaries of the ST64B DX64 prophage in each genome. The 
ST64B sequences in each strain were then selected out. The 
integration point for the ST64B DT64 prophage is located at 2 1 042 



which is between SB27 and SB28 around the midpoint of the 
ST64B D T64 genome. The comparison of the DTI 04 sequence was 
therefore more conveniently done against the sequence from DT2 
which was very close to the ST64B DT64 sequence. This was done 
using BLASTn from the NCBI Align program. Maps showing the 
aligned and non-aligned regions were constructed. The gene map 
from ST64B DT64 was used to determine the positions of the genes 
and intergenic regions for the DT2 sequence so that the 
corresponding positions on the DTI 04 sequence could be located. 
The non-aligned regions and areas of relatively low alignment 
were located and were extracted from each of the DT2 and 
DTI 04 sequences. A visual comparison of the ST64B sequences 
for DT2 and DTI 04 was made using MAUVE [20]. 

A similar process was followed for Gifsy-1. Additional whole 
genome sequences for ATCC 14028 (GenBank Accession 
No. NCJ316856) and SL1344 and D23580 (Sanger site) were 
downloaded. The Gifsy- 1 L t2 sequence Accession No. NC_ 
010392.1 was BLASTed against the DTI 04 whole genome 
sequence to locate the 5' boundary and the following 50,000 
base sequence was downloaded. The LT2 and DTI 04 Gifsy-1 
sequences were aligned and a map was constructed to identify the 
aligned and non-aligned regions. Following the discovery of the 
duplication of ST64B genes in Gifsy-lDTi04 sequence this 
sequence was also aligned with both ST64B DX64 and ST64B DX1()4 
sequences. Regions of the Gifsy- 1dtio4 sequence which failed to 
align with any of the prophage sequences were extracted. In the 
same way Gifsy-1 sequences were extracted from the genome 
sequences for DT2, SL1344, D23580 and ATCC 14028. The 
Gifsy-1 sequence from S. enterica subsp enterica Serovar 4, [5], 1 2:1:- 
strain CVM23701 was extracted from contigs ABAO0 10000 16.1 
and ABAO0 1000067.1 from whole genome shotgun sequencing 
project ABAO00000000.1 BLASTn was used to compare the 
various Gifsy-1 sequences, and regions which appeared to be 
unique to SL1344, DT2 and CVM23701 were extracted. (The 
sequences for D23580 and ATCC 14028 were not significantly 
distinct from that of DT2). 

The co-ordinates for the Gifsy-3 prophage in ATCC 14028 
(1284562 and 1335445) were obtained by aligning the LT2 
genome sequence with a sequence from the ATCC 14028 genome 
estimated to contain the prophage knowing its approximate 
location from the genome map [21]. Gifsy-2 from LT2 was 
downloaded using the published co-ordinates [22]. BLASTn was 
used to locate and extract a region in Gifsy-3 which did not appear 
in Gifsy-2 or any of the Gifsy-1 sequences. Primers were chosen for 
the amplification of a portion of this sequence (Table 2). A visual 
comparison of the Gifsy-1 sequences for LT2, DT104, DT2, 
SL1344 and CVM23701 and the Gifsy-3 sequence from ATCC 
14028 was made using MAUVE [20]. 

MLVA Typing 

For routine MLVA typing the procedure of Lindstedt et al [1 8] 
was followed except for minor modifications: the concentration of 
primers for STTR6 and STTR9 was halved, the size marker was 
labelled with ROX and fragment sizing was performed on an ABI 
3 1 30 sequencer. MLVA profiles were expressed as a number string 
of allele sizes (as generated by the sequencer) for each VNTR in the 
order STTR9-STTR5-STTR6-STTR10pl-STTR3. Additional S. 
Typhimurium VNTRs were amplified using primers shown in 
Table 2. PCR procedure and product detection were as described 
below for CRISPR typing except that PCR conditions of 30 cycles 
of 94°C for 30 s, 57°C to 58°C for 30 s and 72°C for 60 s to 90 s 
with a final 72°C for 7 min were applied and for STTR1, STTR7, 
STTR8, STTR11 and STTR12 5% DMSO was included in the 
mastermix to eliminate non-specific product formation. 
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CRISPR Typing 

The two CRISPR sequences in each isolate were amplified 
using a number of primer pairs depending on the size and spacer 
composition of the CRISPR sequence present (Table 2). For alleles 
up to 1500 bp adequate amplification was usually obtained using 
primer pairs spanning the full CRISPR but for longer sequences 
two PGR products were generated using internal primers derived 
from spacer sequences. The mastermix contained 2 mM MgCl 2 , 
5 pmol of each primer (Geneworks, Adelaide, South Australia) and 
0.5 U of AmpliTaq Gold (Applied Biosystems, Foster City, Calif); 
the initial cycling step was 95°C for 10 min followed by 30 cycles of 
94°C for 15 s, 55°C for 30 s and 68°C for 2.5 m with a final 68°C 
for 7 m ; a 6 aliquot from each PCR tube was electrophoresed in 
a 1.5% agarose gel containing 0.5 ug/ml ethidium bromide at 
80 V for 60 min. PCR products were sequenced using the Applied 
Biosystems Big Dye Terminator v3. 1 and Applied Biosystems 3130 
sequencer. The sequences were analysed by ChromasPro. The 
CRISPR finder program (http://crispr.u-psud.fr/) was used to 
locate the regular repeats and the intervening spacer sequences. 
Results were represented as filled rectangular blocks for 'spacer 
present' or an X for 'spacer absent' in the same order as for S. 
Typhimurium spacers in Fabre et al [23]. 

Amplification of Phage Loci 

For amplification of SB6, SB26, SB28, SB37 and SB46 loci in 
ST64BrjT64 prophage the primer pairs used were those of Ross 
and Heuzenroeder [8]. Table 2 shows primers pairs chosen for 
amplification of seven loci in ST64Bryriu4 which had low or no 
identity with the ST64B sequence from DT2, of five loci in the 
Gifsy-lDTi04 with no identity with either Gifsy-1 L T2 or ST64B 
sequence, of six loci in Gifsy- Isl 1344 an d three loci in Gifsy-lDT2 
with no identity in other Gifsy-1 sequences, of one locus in Gifsy- 
1lt2 shared by Gifsy- Isli 344 an d Gifsy-lDT2 but not by Gifsy- 
1dtio4 or Gifsy-2, three loci in Gifsy-1 cvM2370ij two with no 
identity in other Gifsy- 1 sequences and one shared with all Gifsy-2 
sequences and with Gifsy- 1lt2> an d one locus in Gifsy-3 1402s- 
There are also primers for a complex 22/39 bp tandem repeat 
which occurs in Gifsy-2 but also in Gifsy- lcvM2370i> where the 
repeat sequence is longer. The mastermix composition was the 
same as for CRISPR typing. For amplification the annealing 
temperature was between 55 and 60°C for 30 s and extension took 
place at 72°C for 30 to 60 s. Products were visualised on gel and 
representative ones were sequenced as required. 

Phage Typing 

Isolates were sent to the Microbiological Diagnostics Unit, 
University of Melbourne, Australia for phage typing by the 
Anderson scheme and antibiotic sensitivity testing. 

Supporting Information 

Figure SI VNTR/CRISPR 1 genotypes for panel of well 
characterised S. Typhimurium strains ordered into Major Groups 
and Repeats Groups (RG) based on VNTR (MLVA allele string 
STTR9-STTR5-STTR6-STTR10pl-STTR3+allele sizes for six 
other VNTR) and CRISPR 1 typing (solid block = spacer present, 
X = spacer absent). The spacer positions are same as those applied 
by Fabre et al [34]. The 74 bp spacer includes Spacer 10. Spacer 
22 contains a 6 bp tandem repeat Green = 3 repeats, Black = 4 
repeats, Red = 5 repeats, Blue = 6 repeats. 
(XLS) 



Figure S2 VNTR/CRISPR 2 genotypes for panel of well 
characterised S. Typhimurium strains ordered into Major Groups 
and Repeats Groups (RG) based on VNTR (MLVA allele string 
STTR9-STTR5-STTR6-STTR10pl-STTR3+allele sizes for six 
other VNTR) and CRISPR 2 typing (solid block = spacer present, 
X = spacer absent. The orange block indicates a SNP variant of 
spacer 10 as per Table 6 in [34]). The spacer positions are same as 
those applied by Fabre et al [34] . 
(XLS) 

Figure S3 Diagrammatic comparison (not to scale) of the ST64B 
DNA sequences of S. Typhimurium strains DT2 (Sanger) and 
DTI 04 (NCTC 13348) showing regions of comparative identity 
interspersed with regions of non-identity. 
(XLS) 

Text SI The DT104-specific SB26 sequence in earlier isolates of 
DT197. 

(DOC) 

Text S2 Variation in the ST64B D T64 SB46 sequence for the 

three phage types in RG14. 

(DOC) 

Text S3 Possible origin of the SB46 variant found in the isolates 
in RG6B. 

(DOC) 

Text S4 Alignment of Gifsy-1 DT104 with ST64B DXfi4 and 

ST64B DT104 . 

(DOC) 

Text S5 Alignment of Gifsy- 1 sequences. 
(DOC) 

Text S6 The tandem repeat for SB42 in isolate 06P23331584. 
(DOC) 

Text S7 Variation in the 24 bp VNTR in Gifsy- 1dt2 among 

Repeats Groups. 

(DOC) 

Text S8 Supporting evidence for the identification of ST64B 

sequences in Gifsy-1 D T64 an d Gifsy-1 DT104 . 

(DOC) 

Text S9 The separation of Major Groups 1 and 2 based on 

VNTR data. 

(DOC) 
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